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RADIATION AND COUPLING BETWEEN TWO COLLINEAR 


OPEN ENDED WAVEGUIDES 

By Robert J. Mailloux 
Electronics Research Center 

SUMMARY 

This report describes an experimental and theoretical 
study of collinear coupled waveguides radiating through a 
common perfectly conducting ground plane. The problem is 
formulated as a set of simultaneous integral equations and 
solved approximately by expanding the aperture field in a 
Fourier series. 

The special case of a single isolated waveguide is treated 
first, and the results are compared with experiment as well as 
with the approximate results of Lewin. Attention is given to 
the edge singularities. 

The more general case of two collinear coupled waveguide 
slots is solved by the same numerical method using the 
symmetrical properties of the geometry. Phase and amplitude 
of the mutual coupling are evaluated and compared with 
experiment. Agreement is excellent. 

Considerable attention is given to comparing more simple 
approximations with this result, and on several occasions it 
is pointed out how these can be used and what error the addi- 
tional approximations introduce. 



I . INTRODUCTION 

The radiating waveguide is a fundamental electromagnetic 
structure, and one about which a great deal is known. With 
the realization of large scale microwave arrays, the subject 
of waveguide radiation and mutual coupling has aroused 
renewed interest. Since most of the recent theoretical 
efforts have been concerned with infinite arrays 1, 2,1 and 
since these results are not readily applicable to a considera- 
tion of small arrays, or edge effects in large arrays, there 
is a need to study this problem from the point of view of 
individually coupled elements. This approach has been found 
roost useful in the case of dipole arrays because the results 
can be readily compared with experiment and because the 
analysis may be extended quite naturally to consider the 
edge behavior of small or large arrays and can be used for 
infinite arrays as well. 

The present work is concerned with the study of two open 
ended waveguides which radiate through a common perfectly 
conducting ground plane. The waveguides are arranged 
collinearly as shown in Figure 1, and the structure has been 
analyzed by expanding the solution in a truncated set of wave- 
guide modes as an approximation to the integral equation 
solution. This method was chosen in preference to the compu- 
tationally simpler variational solutions in order that bounds 
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on the accuracy of the theory could be established by 
varying the number of terms used in the solution. The theoreti- 
cal results have been correlated with experimental data in 
order to justify the basic assumptions of the analysis. 

Though many authors have considered the problems of wave- 
guide radiation and mutual coupling , there are still a great 
many questions which have yet to be answered concerning these 
subjects. The following discussion is intended to point out 
some of the gaps in present theoretical knowledge. 

Much of the earlier work in this area was concerned with 
computing the field patterns of single waveguides with assumed 
aperture distributions. The next subject which attracted 
general interest , the input admittance of waveguide radiators, 
has traditionally been studied ' using variational formulas 
with the lowest order mode as an aperture distribution. 

These formulations have been quite successful and give 
reasonably good approximations to both the conductance and 
susceptance. As mentioned earlier, they cannot easily be 
made to yield increasingly more accurate solutions. Other 
stationary formulations^' 7 *8,9 have been studied for cases 
involving waveguides radiating into complex media. These 
formulations are based upon the use of Fourier integral 
techniques and use the dominant mode field as the aperture field 
approximation . 

Galejs 10 was apparently the first to mention that the 
electric field distribution in the waveguide E-plane was not 
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uniform. He correctly observed that the electric field 
normal to the waveguide edges should have a singularity of the 
same order as the static field singularity (r -1 / 1 ) at a 90° 
edge. His solution, which was formulated by considering the 
junction between two waveguides as the size of one was made 
to approach infinity, assumed a trial field which exhibited 
this singularity. He found differences in admittance of 
about 5% when comparing this field approximation with the 
uniform field approximation. 

In contrast to the single waveguide radiation problem, 
relatively little work has been done to describe the coupling 
of two waveguides. It appears that the first study of this 
sort was performed by G. W. Wheeler 11 who assumed the coupled 
radiators were in the far field of one another. By comparing 
theory and experiment, Wheeler was able to show that a single 
mode solution was quite adequate for the waveguides which he 
considered when coupled under this far field condition. In 
1956 Levis 12 derived general equations for a variational 
formulation to obtain the coupling between a number of 
generally cylindrical waveguides radiating through a common 
ground plane. He applied the method to a set of coupled 
annular slots. Galejs 12 applied a stationary formulation due 
to Richmond 14 to solve the problem of two parallel slots in 
a ground plane, with both slots backed by waveguides. He 
also considered the limit in which the slots were the same 
size as the waveguides. His method yielded usable and 
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convenient formulas; it includes the implicit assumption 
that the tangential magnetic field at the coupled waveguide 
aperture is the same as the magnetic field which would be 
present on the ground plane if the coupled aperture were not 
present. In this manner Galejs avoided the problem of 
solving an integral equation. 

The most recent work in this area is described in a 
report by Lyon 15 et.al., and concerns the power coupling 
between various radiating structures including two open ended 
waveguides in a common ground plane, with arbitrary orienta- 
tion and spacing. This experimental and theoretical study 
resulted in much useful data and some very convenient approxi- 
mate formulas were based upon a single mode approximation to 
the coupling and also included the assumption, equivalent to 
that implicit in the work of Galejs, that the waveguide-backed 
ground plane slot has a total magnetic field equal to the 
incident magnetic field. This assumption again eliminated the 
need to solve an integral equation to determine the coupling. 

The recent work on infinite arrays has shown an increased 
awareness of the importance of higher order modes in scanned 
array situations. Galindo and Wu 1 ' 2 have studied a limited 
but important class of infinite arrays for which the boundary 
value problems could be rigorously solved in terms of scalar 
functions. Farrell and Kuhn 3 have used two sets of LSE modes 
to show the existence of a deep null in the array power 
pattern of a triangular grid array as the array was scanned 
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only a few degress off broadside in the H-plane. This null 
had been detected experimentally at an earlier date, but was 
not revealed by theoretical results based upon a single mode 
approximation. The authors also mention an unpublished 
study in which the excitation of cross polarized modes is 
studied as the array is scanned in the E-plane. This appears 
to be the first mention of this phenomenon. 

There are several other mathematical techniques which 
have been found useful for related problems but do not seem 
to have been used for rectangular waveguides. The geometrical 
theory of diffraction has been applied to study the radiation 
of horns^®'^ and the mutual coupling^® of parallel plate 
waveguides. The Weiner-Hopf procedure^ has been used to 
study radiation of circular waveguides. 

This brief survey emphasizes the fact that there are 
many areas which should be studied further in order to 
achieve a better understanding of this basic problem. The 
present work appears to be the first example of a mutual 
coupling study (excluding the two infinite array studies) 
in which an attempt has been made to solve the integral 
equations which govern the coupling of two waveguides. It 
also appears to be the first example in which more than a 
single mode is used to describe the coupling, and in addition 
it presents the only experimental data currently available 
showing both the amplitude and phase of the signal coupled 
between two collinear open ended waveguides. Though this 
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study is intended as an aid in evaluating simpler and more 
approximate analyses , it cannot be used to estimate the cross 
polarized field components excited at the array face. These 
components are not large for the collinear coupled slot case 
discussed in this report and the analysis presented here has 
been found to be an excellent approximation for all cases 
tested experimentally. 

II. FORMULATION OF THE BASIC PROBLEM 

The basic geometry for this study is shown in Figure 1. 

Two rectangular open-ended waveguides are mounted flush with 
an infinite, perfectly-conducting ground plane. The free 
space field in the half space bounded by the perfectly- 
conducting plane (with apertures) can be written in the conven- 
tional way in terms of the tangential aperture field as: 

(for z ^ 0) 

B(r) = j 2 <jj ut E I r° (r , r ’ ) • (z x E)d S' 

p=l,2 S p 

(1) 

E (r) =2 E / VG (r ,r ' ) x (z x E)d S' 

P-1.2 S p P 

where r°(r,r') = (U + A-vvJG (r ,r ' ) 

k o 

e -jk 0 |r-r'| 

and G (r ,r ' ) = , , 

4 tt | r-r ' | 

and | r— r ' | = /(x-x) 2 + (y-y') 2 + z 2 


8 



An exp(+jwt) time dependence is assumed and has been 
suppressed. Vectors are denoted by a bar above the expression 
and dyadics by a bar below. U is the unity dyadic and r° is 
the conventional free space dyadic Green's function. An 
expression which is entirely equivalent to (1) , but which 
helps to explain the major approximation of the analysis to 
follow, is obtained by defining the magnetic hertzian poten- 
tial 


- e - jlc ° I r_r ' 

(z * E) — 


( 2 ) 


T m 


(r) = 


2 TT ( 


p=l,2 S. 


-r- ' 


r-r 


d S£ 


The corresponding fields are written 

(3) B(r) = v(v • 7 m ) + k£ 7 m 

E (?) = -j b) V x n m . 

The above equations define the fields in free space bounded 
by the conducting plane with apertures. 

The field in each of the rectangular waveguides may be 
written in a number of ways, but for the purposes of this 
development it is convenient to expand the field in terms of 
the magnetic hertzian potential as done for the free space 
region. This is done by defining two scalar hertzian poten- 
tials iTmx and w m y such that 

(4) 7 m = x Tt mx + y 7r my . 
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Expressing the fields in this form is convenient because 

it leads to a very simple and symmetrical form for the 

integral equations, and also because it emphasizes the formal 

similarities between this problem and several well-known 

(20 21 22 23 ) 

aperture diffraction problems tit . This expansion 
is clearly valid for the exterior half space, and one might 
expect intuitively that it would also be valid within each of 
the waveguides, where it is equivalent to two sets of LSE 
modes. These sets are not sufficiently complete to expand 
the solution of an obstacle in the waveguide, but it is 
proven in the appendix that they are adequate to expand the 
waveguide fields when the matching takes place at a plane 
perpendicular to the waveguide axis. If the waveguide aper- 
ture field is used in equation (2) , it is only necessary to 
equate the tangential magnetic fields on both sides of each 
aperture in order to obtain integro-dif f erential equations for 
the coupled fields. In order to distinguish between the free 
space fields and the waveguide fields the notation tt ' , n ' 

will be used to designate the free space fields, while the 
un-primed potentials will be used (later with subscripts) to 
designate the fields in each waveguide. With these definitions, 
the equations that express the continuity of the magnetic 
fields in the apertures can be put into a form similar to that 
used in the diffraction problems mentioned earlier. At each 


aperture 


( 4 ) 


+ 


9x z 9y 


- 3 - + k 2 ) u; 


mx 


1T mx ) 


= 0 


( 5 ) 


(_li + 

9x 2 9y 


- + k ) (n ’ 

2 o my 


TT ) 

my 


= 0 


( 6 ) 


9x 11 my) 


9y 


^mx "mx^ 


These equations have been solved rigorously in only a very 
small number of cases, most of which are mentioned by Bouwkamj^®) . 
The radiating waveguide is not among these. All of the 
published results describe some approximate solution to the 
above set of integro-dif f erential equations. The analysis 
which follows is also based upon a simplification of these 
equations, to the extent that a single hertzian potential 
function will be used to approximate the input admittance of a 
single radiating waveguide and the coupling parameters of two 
collinear radiating waveguides. This assumption is not new; 

it is implicit in the work of Lewin and Edelberg and 

( 24 ) ( 23 ) 

Oliner' and it is made explicit by Kieburtz and Ishimaru 

The study presented here differs from others studies making 
the same basic assumption in that no other major assumptions 
are made. Instead, the reduced integro-dif f erential equations 
are solved by expanding the solution in a truncated infinite 
series of a restricted set of modes, and then inverting the 
matrix equation which results. The convergence of this tech- 
nique has been studied quite carefully and the truncation 
errors are estimated whenever the data is quoted. By these 
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means, this theoretical study results in a very accurate 
solution of the approximate model used. Attention is then 
focused upon the comparison of this model with the experimen- 
tal results to verify the validity of the model at hand. 

The assumption that the single hertzian potential function 
7T mx suffices to expand the field ultimately requires that there 
be no component of electric field in the x direction. Equations 
(4) and (6) reduce to the following: 


171 ‘lb- + k o> ( '« - " 0 


i7 ( ’mx ‘ "mx* = 0 


At each aperture a solution of these equations is 


(9) *mx = "mx + A cos k o x + B sin k o x 

(4) 

This integral equation was first presented by Lewm 

without mention of its approximate nature. The approximation 

can be understood by noting that the factors A and B are 

constants and not functions of the spacial parameter y. 

( 21 ) 

Copson has pointed out that the constants of integration 

must be determined so that the fields tangential to an edge 
vanish at the edge. This is necessary to assure that all 
fields are square integrable in a three dimensional domain 
around the edge. The constants of integration for this 
problem must therefore be determined so that the tangential 
alectric field vanish at the edges x = + a / 2 * This condition 
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is the Babinet equivalent of that ordinarily encountered in 
dipole theory (25) states that the current must vanish 

at both ends of the dipole. However, though this condition 
is satisfied by a constant for the infinitely thin dipole, 
the two constants available in this solution are only capable 
of assuring that this condition is satisfied at one point at 
each edge. Therefore it is clear that the solution used in 
this analysis is rigorous for slots whose width "b" equals 
zero, but that it is only approximate for the present case 
with finite slot width. The edge condition is not violated 
in this study however, because the solution is expanded in 
terms of a finite number of waveguide modes, each of which 
satisfies the stated condition at x = + / 


III. SOLUTION FOR THE SINGLE RADIATING WAVEGUIDE 
Subject to this approximation, the electromagnetic field 
is written below in component form. 

( 10 ) 


3 2 7T 


E x = 0 


B x = 


mx 


3x 2 + k o "mx 


_ . 8 "mx 

E y ‘ 


3 2 TT 

B = 

y 3x3y 


r _ .. 3lI mx 

E = + ju> 

z ay 


3 "mx 
B z 3x 3 z 


The fields in each waveguide are obtained by writing the 


total potential function in terms of an incident mode field 
plus an infinite number of reflected waves. The incident 
component of the hertzian potential is normalized to unity. 
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( 11 ) 


mx 


y z 

= (e 10 - 


- y Z X+ a / 2 

re 10 ) sin 7T ( — ) 


+ Z 
m=l 


OO | 

^ A mn 
n=0 


sin nnr ( 


x + a /2 } 

a 


cos nn ( ^* — ) 


y z 
e mn 


(12) where Ymn = + /[ (2f) 2 + (“)2 ~J -~k£ 

( <* mn for (m,n) ? (1,0) 
j^-je iQ for (m,n) = (1,0) 


The prime symbol above the summation is used to indicate the 
omission of the (m,n) = (1,0) term from this double summation. 
The y component of electric field is given by (10) and (12) 
and given below at z — 0 


3 7T 


y 1 

y z=0 


= -]W 


mx 


3 z 


z=0 


= -jsi J 


V 1 


o 


(1+D sin tt (• 


x+ a /. 


') 


(13) 


• x+ a / 2 y +b /2 

* E „ A mn ^n sin m7T > cos n " ( ~ b“ 1 1 

m=l n=0 


This expression is no longer a convergent fourier series 

because the electric field perpendicular to the top wall 

-1/3 

of the waveguide exhibits the p singularity typical of a 

90° edge, where p is the distance from the edge. The 
expression (11) , however, is convergent everywhere in and out- 
side of the aperture, and this property is required to assure 
the validity of this solution. Singularities of this same 
order are exhibited by E z at y = + b /2 and by B z at x =t. a /2 • 
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The basic equation (9) is re-written below using the 
given hertzian potential form inside of the waveguide. 

(14) 



a/2 

/ 

-a/2 


b /2 e -jk o r 

f E ( x ' , y ' ) — — dx 1 dy ' 

^ /o y r 


A cos k Q x + B sin k Q x 


X+^/ n oo 

+ ( 1— r ) sin tt ( — ) + E 

a m=l 


i co x+ a /? y + b /2 

E A mn sm mu ( — ) cos n tt ( , ) 

n=o a D 


The solution proceeds by using the electric field from 
(13) in the left side of (14) , and setting B = 0 because of 
the even symmetry. After defining the integrals: 

(15) 

1 a/2 b/2 e~ jk o r n x '+ a / 2 y' + V2 

I (x,y) = — / / sin mu ( ) cos mr ( H ) dx 1 dy 1 

2 " -a /2 -b /2 r n 


where 


: l l = ^ ^ x-x ' 


) 2 + (y-y' ) 


one obtains the reduced form of equation (9) applicable to the 
single radiating waveguide. 


l o 


(r+l) I 10 (x,y) - 


? 


m=l , 3 n= 0 , 2 


I 0 A mn Y mn I mn ( X *Y) 


(16) 


= A cos k Q x + (1-r) sin n ( X+ J™ -) 


+ E I A sin imr ( x+ Zj. ) cos nn (^— 7 -——) 

m=l , 3 n= 0 , 2 mn a b 

The set of modes A mn is truncated at finite values of M 
and N (M T ,N T ) , and the constant A is treated as an unknown. 
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This determines a matrix equation when the equation is satis- 
fied at a number of points within the aperture. The number 
of points chosen is one greater than the number of modes 
used in the solution. Treating the constant in this manner 

is equivalent to a procedure which is commonly used in antenna 
(25) 

theory , namely that of defining a difference potential and 
assuring that it is zero at some point on the antenna. This 
method is applicable when the chosen distribution functions 
of current (or in this case electric field) are such that the 
fields automatically satisfy the edge condition. One row of the 
matrix equation is shown below. 


t-Y 


1 0 


1 0 


x + a /2 

(x,y) + sin u (- -) ] r-A cos k x 

d u 


m t *n t 


(17) - Z Z [y 

m=l,3 n=0 , 2 


mn 


I mn (x ^ ) 


,x+ a / 2 . 
+ sin nm ( ) 

a. 


y+ b /2 

cos nTT (■* — ) ] A 

b mn 


x+ a / 2 

= y 10 I 10 (x , y ) + sin it ( ) 

a 

This equation was repeated at points on a rectangular 
lattice within the aperture. 


IV. SOLUTION FOR TWO COUPLED COLLINEAR ANTENNAS 
When two collinear slots couple, the coupled integral 
equations can be reduced to two uncoupled integral equations 
by virtue of the geometrical symmetry. The basic integral 
equations are written below. 
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At Aperture #1 


a/2 a/2 e _jk o r i l 

/ / E v (x',y') - dx’ dy’ 

y l x r, x 


2ttoj -a/2 -b/2 


l 1 


(18) 


. a/2 b/2 -jk Q r 12 

7- 1 / 1 E y 0 (x 2' y,) — 7 dx 2 dy ’ 

2tIuj -a/2 -b/2 2 12 


A cos (k 0 x^) 


x,+ a /2 

+ B sin k Q x^ + S 1 (l+r 1 ) sin tt ( — i ) 


,1 


00 00 1 . X1+V9 v+^/o 

+ £ £ Si A__ sin mw (— ± — ) cos mr — -^- 2 -) 

ra=l n=0 a b 


At Aperture #2 


. a/2 b/2 „-jkr 2 , 

/ / E y (x| ,y ' ) e <2 21 dx’ dy’ 

2 ™ -a/2 -b/2 Y1 r 2 i 1 


a/2 b/2 


2lTtu -a/2 -b/2 y 2 2 


-3 k o r 22 


2 2 


dx 2 dy' = C cos(k Q X 2 ) 


(19) 


^ / 2 

+D sin k Q x 2 + S 2 (l+r 2 ) sin n ( — ) 


+ £ 


, x ? +a /2 


00 00 A n ’ / y t 

S 9 B 

m=l n=0 mn 


s 2 b mn sin mu (— ^ — — -) cos m ( y+ t./ 2 ) 


where r u = / (x^x^) 2 + (y-y') 2 r )2 


r 22 = / (x 2 -x 2 )? + (y * y ' )J r 2 1 


/ (x^+D-x^ 2 + (y-y 1 ) 2 

/ (x 2 +D-x|) 2 + (y-y 1 ) 2 
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These equations can be uncoupled by considering separate 
symmetric and anti-symmetric excitations S 2 = +S^ = +1. When 
this is done there is complete symmetry or assymmetry about 
the line x 1 = d /2- Setting -x 2 = +Xj^ = x, C = +A, D = TO and 

B mn = ± A mn for m odd and B mn = '* rA mn' for m even ' one can see 
that the resulting equations have the solution Ey 2 (-x 2 ,y) = 

+ E (x-, ,y) and that the two equations reduce to a single 
“ Y 1 

equation in these two limits. It is then a simple matter to 
reconstruct the case of arbitrary excitation by taking the 
appropriate combinations of these two cases. Later the 
solution to the case of parasitic antenna #2 (S 2 = 0) will be 
considered. 

Defining 

(20) A 

I mn < x 'y> = I mn (x ' y) ± I mn (_x + d ' y) 
and using the previous definition (15) of this integral, one 

obtains the following equations for the symmetric and anti- 
symmetric problems. 

S , s 
A x+ /2 A 

[-Y 1 o I 10 ( X 'Y) + s i n 77 ( — :j— : -) ] r - A cos k Q x - D sin k Q x 


( 21 ) 

S b S 

oo 0 ° J\ x+ a /? y+ D / 9 A 

- Z Z [ y I _ (x,y) + sin mir ( -) cos nn ( £ K ) ] A^ n 

m=l , 2 n=0 , 2 mn mn a b "mn 

S 

■r A , x . • ,x+ a /2x 

= Y io I io (x ' y) + sin 7,( — 1 
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As before, this equation was written at a finite number 
of points within the aperture and a truncated series of modes 
was used to describe the field. There are now two constants 
per equation, and these are determined as before by requiring 
the integral equation to be valid at two extra points in the 
aperture . 

In order to construct the solution for the case when 
antenna #1 is excited and antenna #2 is kept parasitic and 
terminated in a matched load, it is simply necessary to super- 
impose the symmetric and antisymmetric solutions and to 
divide by two in order to normalize to the source amplitudes. 
Following the usual procedures for such problems and using 
the subscript p to denote this case, one obtains the 
coefficients below: 

( 22 ) s x j = s 22 - %(r s + r A ) 

S A 

(23) s 12 = s 21 = Mr - r ) 


(24) A mn 
P 


(25) 



Hi A, 


mn 



+ A 


mn 


- B 


mn 


) 

) 


V. THEORETICAL AND EXPERIMENTAL RESULTS 
The convergence of the matrix solution to equation (21) 
has been studied in detail in order to determine its accuracy. 
Furthermore, since this entire mathematical development rests 
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upon the applicability of the scalar approx imat ion (9) , an 
experimental program was undertaken to provide results to 
compare with the theory. These two subjects are discussed in 
detail below. 

The convergence properties of the solution for single 
and coupled waveguides are basically different and so they 
are discussed separately. Table I shows the variation in 
input reflection coefficient for the single radiating wave- 
guide as various sets of modes are chosen to represent the 
field distribution. All accuracies are relative to the 
highest order solution (Case 6) . The single mode approxima- 
tion yields an imaginary part which is in error by about 
twenty percent and a real part (much smaller than the 
imaginary part) which has the wrong sign. A solution with 
four modes (Case 3) but including the 12 and 14 modes to aid 
in approximating the E-plane distribution is among the more 
accurate solutions shown. This high accuracy is coincidental 
and was not generally repeated, but this solution was always 
within about ten percent of the highest order solution for 
all cases examined, and so was later used to form the basis 
of the coupled waveguide study (after the addition of another 
mode having anti-symmetrical x-dependence ) . Another four 
mode solution (Case 2) with one less mode to describe the L- 
plane distribution was less accurate than the results of 
Case 3. Nine modes (Case 4) serve to determine the imaginary 
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TABLE I 


A. SINGLE WAVEGUIDE REFLECTION COEFFICIENT 

a/ A = 1.0; b/A = 0.3 

CASE MODES REFLECTION COEFFICIENT 


1 


10 




r 1Q ■ -0.3928 - 

jO. 32680 

2 


10, 12, 30, 

32 



r 10 = 0.10230 - 

jO. 24786 

3 


10, 12, 14, 

30 



r 1Q = 0.07742 - 

jO. 26969 

4 


10, 30, 50, 

12, 

32, 


r 10 = 0.08309 - 

jO. 27193 



52, 14, 34, 

54 




5 


10, 30, 50, 

70, 

12, 

32, 

r 10 = 0.08319 - 

jO. 27085 



52, 72, 14, 

34, 

54, 

74 


6 


10, 30, 50, 

12, 

32, 

52, 

r 10 = 0.07729 - 

jO. 27786 



14, 34, 54, 

16, 

36, 

56 




B. MODE 

AMPLITUDES 

FOR CASE 6 


r io 

= 

0.07729 - jO. 27786 




A ^ “ —0.01965 

- jO. 03589 

A 30 

s 

-0.05521 + jO. 00195 




A 34 = -0.00302 

+ jO. 00051 

a 50 

= 

-0.00638 - jO. 00151 




A 54 = -0.00076 

- jo. 00034 

a 12 

= 

-0.05538 - jO. 08904 




A 16 = -0.00736 

- jO. 01418 

a 32 

= 

-0.00762 + jO. 00171 




A 3 g = —0.00107 

+ jO. 00030 

A 52 

s 

-0.00165 + jO. 00022 




A5g = —0.00021 

+ jO. 00021 



part of the reflection coefficient to within about two percent 
of the highest order solution considered and the real part to 
about eight percent. A comparison of this solution (Case 5) 
with Case 6 illustrates that the distribution in the Il-plane 
(x-direction) is fairly well approximated by a few terms 
because the three extra modes introduced to more precisely 
approximate the x dependence barely change the value of the 
reflection coefficient. Table IB lists the mode amplitudes 
for another twelve mode solution (Case 6) which was chosen to 
provide the highest order approximation to the E-plane distri- 
bution, and hence the most accurate solution for the 
reflection coefficient. Inspection of these mode amplitudes 
is a convenient measure of convergence because the integral 
equation (21) is solved at a number of discrete points within 
the aperture and therefore one can expect an error in the 
reflection coefficient which is approximately the size of the 
modes truncated. Table IB therefore shows why it is especially 
necessary to include the modes which approximate the E-plane 
field variation (12,14 etc.). 

Figure 2 shows the aperture field variation provided by 
this twelve mode solution. The approximation of the singular- 
ity at y = +° / 2 is evident in this distribution, as is the 
fact that the H-plane distribution is fairly accurately 
represented by a single mode. 
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Normalized Electric Field 
Amplitude |Ey| 



FIGURE 2 APERTURE ELECTRIC FIELD 
DISTRIBUTION FOR RADIATING WAVEGUIDE 
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Figure 3 displays a set of curves of input admittance of 
a radiating rectangular waveguide. These curves are presented 
in the form of the well known results of Lewin^ and a few 
points are shown to illustrate the correlation with Lewin's 
data. This correlation is very good except for deviations 
of about ten percent in the conductance for thick slots, and 
of about thirteen percent in the susceptance for very thin 
slots. There is no experimental evidence to compare with 
these two theoretical results, but since they are both derived 
using the same scalar approximation (9) and since the present 
analysis is a convergent solution, it is claimed that the 
data of Figure 3 is a more accurate solution to the integral 
equation. Whether or not these curves show closer agreement 
with experiment will determine the usefulness of the scalar 
approximation itself. 

Regarding the usefulness of this scalar approximation for 

describing waveguide admittance, Figure 4 compares the 

(26,27) 

experimental data of Venema with the results of this 

theory and the work of Lewin. These results indicate that 
the solution in terms of this single set of hertzian potentials 
is indeed a valid approximation for describing waveguide 
admittance. Venema ' s data does appear to follow the results 
of this theory more closely than Lewin's theoretical results, 
but for the size slots chosen there is so little difference 
between the two that no conclusion can be drawn except that 
they both agree very closely with experiment. 
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It is appropriate to point out that Lewin's stationary 
solution, which uses a single mode, is considerably more 
accurate than the single mode approximation used with this non- 
stationary solution. The advantage of the method presented in 
this note is that enough higher order modes can be conveniently 
included to bring about a convergent solution. 




FIGURE 3 ADMITTANCE OF SINGLE 
RADIATING WAVEGUIDE 
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Normalized Susceptance And Conductance B/yo, G/y 0 



a /\ 

FIGURE 4 COMPARISON BETWEEN THEORY 
AND EXPERIMENT 
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The solution of equation 21 for the case of two collinear 
coupled waveguides involves some of the same convergence 
properties as the single waveguide, but it also offers the 
possibility of relaxing the convergence requirements for the 
coupled parameters. The self admittance terms necessarily 
behave in a very similar manner to the single waveguide, and 
so all of the comments about truncation apply here as well. 

The presence of the second slot in the ground plane (with the 
second waveguide terminated in a matched load) causes a small 
change in the input reflection coefficient of the driven 
antenna. Figure 5 shows that the imaginary part of the reflec- 
tion coefficient undergoes a change of about eight percent 
with spacing variations, while the real part being much 
smaller, changes by about fifty percent. The effect of this 
coupling is no longer visible when the center to center dis- 
tance exceeds about 1.2A. When the spacing ^/ \ exceeds five, 
the subtraction indicated in equation (23) begins to intro- 
duce numerical inaccuracies , but for the range of spacings 
used throughout this report no additional approximation is 
introduced by this symme tric-antisymmetric formulation. 

A more familiar property of closely spaced radiating slots 
is the direct coupling between the radiating signal and the 
parasitic antenna. Table IIA shows that for the case considered 
this parameter is relatively insensitive to the approximation 
used in computing the aperture field. 
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Reflection Coefficient Sn 



FIGURE 5 REFLECTION COEFFICIENT OF DRIVEN 
ANTENNA WITH PARASITIC ANTENNA MATCHED 
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TABLE II 


A. COUPLING COEFFICIENT FOR COLLINEAR COUPLED WAVEGUIDES 


a = 2.286 cm 





d 

= 2.921 


b = 1.016 cm 





frequency = 

10.0 GHz 

CASE 


MODES 



COUPLING COEFFICIENT S 12 

1 

10 




s 12 

= 

0.01669 + 

jO. 00754 

2 

10, 

20, 30, 40 

, 12, 

14 

s 12 

= 

0.01819 + 

jO. 01124 

3 

10, 

20, 30, 12 

, 14 


S 12 

= 

0.01786 + 

jO. 01120 

4 

10, 

20, 30, 12 

, 22, 

32, 

s 12 


0.01874 + 

jO. 01047 


14, 

24, 34 









B . MODE 

AMPLITUDES 

FOR CASE 4 


S 11 * °* 

05595 - 

jO. 22203 



s 12 

= 

0.01874 + 

jO. 01047 

A 20 = °* 

00233 - 

jO. 00210 



b 20 

ss 

-0.00294 

- jO. 03949 

P 





P 




0 

1 

(i 

a 

m 

< 

.02432 

- jO. 00741 



B 30 

= 

-0.00344 

+ jO. 00266 

p 





P 




a 12 = 

.05513 

- jO. 09512 



b 12 


0.00211 - 

jO. 00130 

P 





P 




A 22 = °* 
P 

00030 - 

jO. 00033 




— 

-0.00085 

-jO. 00569 

a 32 = ”° 

.00652 

- jO. 00053 



b 32 


-0.00092 

+ jO. 00095 

P 





P 




a 14 = 

.01430 

- jO. 02800 



b 14 


0.00058 - 

jO. 00041 

P 





P 




A 24 ~ 0 . 

00004 - 

jO. 00010 



b 24 


-0.00025 

+j0. 00161 

P 





P 




A£4 - -0 

.00184 

- jO. 00009 



b 34 

P 

— 

-0.00026 

- jO. 00027 
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Table IIB shows the coupling coefficients for Case 4 with one 
antenna driven and the other terminated in a matched load. 

This data is revealing because it demonstrates that the 
error involved in estimating due to approximating the 

aperture field singularity is likely to be significantly 
larger than that involved in evaluating S 21 for the same field 
approximation. (Notice that the ratio of | I to l^ni is 
much larger than the ratio of | I to l s 12l)* This is 
because the entire field exhibits this singularity, not merely 
the reflected wave. Table IIB does show that the ratio 
j a 12 I/I l +s li I ^- s indeed about the same (0.1) as the ratio 
IB 12 i/1 S i2 1= (0.115) . 

These comments relate to the convergence of the chosen 
approximate solution, but reveal nothing about the correspon- 
dence of this analysis to the physical problem of collinear 
slot coupling. It remains to show that the solution based 
upon this approximate scalar model has relevence to the 
problem studied. 

An experiment has been conducted to test the validity of 
the analysis for the collinear ly coupled waveguide case. 

Since the coupled signal was so small compared to the trans- 
mitted signal, none of the conventional techniques for 
measuring the phase of the coupled signal sufficed. To solve 
this problem a new technique was devised which made it 
possible to measure the phase to within +4°, even when the 
coupled signal amplitude was -40 db relative to the transmitted 
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signal. Figure 6 shows that this data compares very closely 
with two theoretical solutions. The five and nine mode 
theoretical solutions yield nearly the same coupled power, 
but the difference in phase between these two approximations 
can be as much as 15°. This discrepancy is caused by the 
omission of the higher order modes (22) and (24) in the 
five mode solution. These modes possess odd symmetry about 
the center of each waveguide and become very large at the 
high frequencies where the solutions differ substantially. 

The amplitude of &22p a bout one third Is-^l at 10.0 GHz, 
but it is about equal to | S 2 I at 12.5 GHz. 

Since the five mode solution does yield reasonably 
accurate results and requires much less computation time, it 
has been used in Figures 7a and 7b. These figures show the 
amplitude and phase of the coupled signal for waveguides with 
the dimensional ratios a /b = 2.0 and 2.5 and for waveguide 
widths up to a / X = 0.99. The sets of curves with different 
a /b are essentially parallel and sufficiently close together 
to allow interpolation between the. Therefore, the curves 
provided cover the range of almost all commercially available 
wave guides. 

Several important points are noted when these results 
are considered in detail. The first relates to the work of 
Lyon et.al.^"^ in which it is observed that the coupling 
amplitude Is^l varies very nearly like as the distance 

d between slots is varied, and that this is approximately 
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Frequency (Ghz) 

A. Coupled Power - Theory Vs. Experiment 


FIGURE 6 COMPARISON BETWEEN THEORY AND 
EXPERIMENT- COUPLED WAVEGUIDES 
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Phase Of S12 (Degrees) 
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Coupling Amplitude IS12I 



FIGURE 7 COUPLING BETWEEN COLLINEAR SLOTS 
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true even for very close spacing. Unfortunately, this is 
only true in a qualitative way and the coupling does not 
follow this law closely enough to allow its use quantitatively. 
The complex coupling coefficients S-^ does vary like 
[exp(-jk Q d) ]/d 2 for the collinear case when the spacing is 
sufficiently large, but Figure 7a shows a pronounced curvature 
of the amplitude plot, especially for the larger slots 
considered, and so it is obvious that the asymptotic behavior 
of | S 2 1 doe s not begin to dominate unitl x exceeds 2.0 or 
2.5. The convergence of the curves in Figure 7b shows that 
the phase relationship is even less predictable and, though 
this is not shown in the Figure, the asymptotic type of 
variation in phase does not dominate until \ exceeds about 
5. This same phase nonlinearity is exhibited by the more 
accurate data of Figure 6b, and so is clearly not caused by 
the approximations involved in the computations for 7b. 

VI. CONCLUSION AND COMMENTS ABOUT ARRAY THEORY 
This study has shown that a set of fields derivable from 
a single component of the hertzian magnetic potential 
function can provide an extremely accurate solution to the 
boundary value problem for a single waveguide radiating 
through a perfectly conducting ground plane or for two 
collinear radiating slots. The solution is rigorous for 
idealized slots of zero thickness, and is a good approxima- 
tion even when the slot is the size of a standard commercial 
waveguide . 
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This basic formulation could be used to describe the 
mutual coupling between elements of an array, especially if 
the waveguide height is less than one half wavelength, but 
it can yield no information about the cross-coupled fields 
that are excited at the array face. 

The convergence properties of the solution have been 
emphasized, and no attempt has been made to simplify the 
numerical solution. An extensive experimental program was 
undertaken in an attempt to measure the amplitude and phase 
of the coupling as accurately as possible so that these could 
be compared with the theoretical results. The agreement 
between theory and experiment is excellent. 

In addition to the detailed study of the problem at 
hand, this work has led to certain qualitative conclusions 
which have application to the development of an approximate 
array theory for radiating waveguides. These conclusions 
concern the errors which arise using various sorts of 
approximations, which errors seem to fall naturally into two 
categories: those which have to do with specific approxima- 

tions to the aperture fields or the type of solution used 
(ie. a variational theory or the truncated series approxima- 
tion used here) , and those fixed by the physical process of 
coupling independent of the mechanics of solution. Errors of 
the first type have been discussed in sufficient detail in 
Section V, and some hints have also been given in that section 
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relative to simplifying the theory. Concerning the physical 
processes which must be accounted for by any theory, the 
following conclusions now appear evident. 

a) If the waveguide dimensions are such that a cross- 
polarized TE 1q mode (ie. one with its electric field in 
the x-direction) can propagate, then the array theory 
must allow for coupling between the two senses of 
polarization. 

b) This analysis has shown that coupling between collinear 
slots is so small that the reflection coefficient at 

one slot, with the waveguide of the second slot termina- 
ted in a matched load, is very nearly the same as for 
the isolated slot. This is not true for closely spaced 
parallel slots because they couple much more strongly, 
and so it will be necessary in general to account for 
the array geometry in the development of the slot self 
admittance . 

c) The asymptotic approximation S 12 = c[exp (-jk Q d) ]/d 2 
provides convenient insight into the mechanism of 
coupling , but it does not provide accurate amplitude and 
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phase information if the spacing between antennas is 
less than five wavelengths. An accurate waveguide 
array theory must therefore be a near field theory. 
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APPENDIX 


Expansion of the Waveguide Fields in terms of two 
Magnetic Hertzian Potentials. 

The complete field in a rectangular waveguide can be 
expanded rigorously in terms of one component of both the 
electric and magnetic hertzian potentials when these compo- 
nents are both taken along the same axis. The most common 

A 

choice of axis is the direction of propagation (z) but it is 
sometimes convenient to choose these vectors along either 

A A 

the x or the y axis. When the fields are to be matched at a 
plane boundary at the end of a waveguide, the fields can be 
expanded in terms of two magnetic hertzian potentials . To 
show this, it is convenient to choose the magnetic and 
electric hertzian potentials xu^ and xir ex as the basis for 
a rigorous expansion of the waveguide fields for the general 
waveguide discontinuity problem. The proof will then proceed 
by showing that a different set of two magnetic hertzian 
potential functions xir^ and yir m y can be made to satisfy the 
same boundary conditions as the original set. 

An exp[+jwt] time dependence has been suppressed. The 
general electromagnetic field is written below as derived 
from magnetic and electric hertzian potentials. 

(Al) E(r) = -juV - x ^ + v(V • ¥ e ) + k£iT e 
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(A2) 


— k 2 

B(r) - v<v • T m ) + kg.,,, + j-2 V X T e 

As pointed out earlier, the general field inside of a 

A 

rectangular waveguide can be written in terms of xir^ and 
**ex ass 

A ^ 

(A3) E(r) = -.jwV x (xit^) + V(V • xn ex ) + x k 2 Tr ex 

k 2 

(A4) B (r) = V(V • xir mx ) + X k^mx + j-jjj - V x ( XlT 0X^ 


The alternate set of magnetic potentials is given below. 
The prime symbol is used to designate that the magnetic 
potential w ' in this expansion need not be the same as Tr mx 
in the expansion of equations (A3) and (A4) . 


(A5) E(r) = -jwV x (xir^) — j o>v • (y&my) 

(A6) B(r) = V(V • x*^) + V(V • yw my ) + x + y k 2 ir my 


The electric fields tangential to the aperture plane are 
both equated to the tangential aperture field, and therefore 
are equated to each other. In component form, the two 
resulting equations are: 


(A7) 


jw ( 


dir 

mx 
3 z 


3 it 


3 2 IT 


mx 


TF 


) = 


ex 


3x3y 


and 

(A8) 


3 2 Trex 

3x 2 


+ 


+ jw 


3 TT 

my 
3 z 


42 


■ii i mi ■ ■■ ■■ 


II I ■ II II II I I 



If the field can be expanded in terms of two magnetic 


potentials, then the above relations which guarantee the 
equality of the tangential components of £ in the aperture, 
should also guarantee that the tangential components of B 
as expanded using either set of functions is also equal at 


the aperture. 

The magnetic fields tangential to the aperture plane are 
equated to the tangential aperture field and hence to each 


other. The resulting equations are: 


3x3y 


Jr 2 3 u 

3 2 , ° ex 

3x3y *mx * “ 3z 


3 2 tt 

^ + k 2 ir 

3y o my 


(A10) 


+ k 2 ) (ir 
3x 2 ° mx 



3x3y 


If these equalities can be proved to be identically true 
using (A7) and (A8) then the desired proof will be accomplished. 
This can be shown to be the case for equation (AlO) by taking 

2 2 3 

— - — of equation (A8) and by subtracting ju> — of equation 
3x3y 3z 

(AlO) from this. The resulting equation can be written in the 



This equation is identically satisfied using (A7) . 
Similarly, equation (A9) can be shown to be identically satis- 
fied by starting with equation (A8) and using the fact that 
ir ex satisfies the scalar helmholtz equation to re-write (A8) 
in the form shown below. 
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(A12) 



ex 


3ir 


-]u- 


my 


3 z 


Taking + ^-(y^ 2 - + k£) of this equation and subtracting it 

3 

from — of equation (A9 ) , simplifying and again using the 

d Z 

helmholtz equation one obtains the resulting equation: 


3 TT 


(A13) 


3x3y 


[ ju(‘ 


mx 


3Z 


3 Z 


)J = 


a 2 *ex 


3x3y 3x3y 


This equality is also identically satisfied by equation 
(A7) and therefore the two expansions are equivalent for the 
problem considered. 
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